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Abstract. 

We present a simple model for the apparent superluminal 
motion along curved trajectories of features observed in the 
compact radio jet of the quasar 3C 345. The model is inspired 
by the magneto-hydrodynamic approach of Camenzind (1986), 
and assumes a conical geometry of the jet within about three 
milliarcseconds of the radio core. Simultaneous conservation of 



three out of four physical parameters is used to derive ana- 



lytic solutions for the equation of motion of the idealized jet 
j^^ j features. These are the jet's kinetic energy, angular momen- 
l» turn, and momentum along the jet axis. Conservation of an- 
I gular momentum is required to fit the observed trajectories at 
distances larger than 3 mas, in agreement with the kinematic 
predictions of Camenzind's model. The best-fit model variant 
preserves angular momentum, kinetic energy, and opening an- 
gle, with an inclination of the jet axis to the line of sight 
^— ' 8 = 6.8°, an opening angle ij> = 0.95°, and Lorentz factors 



in 
o 
in 



7 = 5.8 and 7 = 4.6 for components C4 and C5, respectively 
(Ho = lOOkms -1 Mpc -1 ). The model ignores turbulent small- 



scale motion of the jet components, and it does not take into 

— «^ account the overall bending of the jet beyond about 4 mas, re- 

w— stricting it's applicability to the inner few milliarcseconds from 

Qnthe core. 
I 

o 

5-H Key words: Galaxies: individual: 3C 345 - Galaxies: jets - 
^3 Galaxies: nuclei - quasars: general - Galaxies: quasars: indi- 
03 vidual: 3C 345 



1. Introduction 

Curved radio structures and bent trajectories of moving fea- 
tures are frequently found in radio jets on parsec and kiloparsec 
scales. In some cases, the bends are quasi-periodic, suggesting 
a helical nature of the underlying structure. Examples are fila- 
ments in M 87 (Biretta et al. 1983), an oscillating ridge line in 
3C 273 (Zensus et al. 1988, Krichbaum & Witzel 1992), bending 
and acceleration in 3C 84 (Krichbaum et al. 1992 and 1993b), 
a bent ridge line in 1803+78 (Krichbaum et al. 1994, Steffen 
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1994) and 3C 380 (Kus 1993). Periodic optical variability, e.g., 
in 3C345 (Schramm et al. 1993) and PKS 0420-014 (Wagner 
et al. 1995), and swings in the polarization angle of OJ287 
(Kikuchi et al. 1988) and 0917+624 (Quirrenbach et al. 1989, 
Qian et al. 1991) have also been attributed to motion in helical 
jets. 

Theoretical models for such oscillating, bent structures in- 
clude helical modes in hydrodynamic jets (Hardee 1987, Owen 
et al. 1989) or in magnetized jets (Konigl & Choudhuri 1985). 
Camenzind (Camenzind 1986, Camenzind & Krockenberger 
1992) developed a model for compact jets with bulk plasma 
acceleration along helical magnetic field lines, based on magne- 
tized accretion disc winds. This model predicts helical motion 
of off-axis features in the jet; through conservation of angular 
momentum in a conical jet geometry the model predicts also 
the asymptotic dampening of the jet oscillations seen in some 
cases. 

The quasar 3C 345 contains such a bent core-jet structure on 
parsec-scales, and also is one of the best-studied sources with 
apparent superluminal motion (Biretta et al. 1986; Zensus et al. 

1995) . Several moving features have been traced during their 
evolution; these 'components' separate from the core on ap- 
parently different curved trajectories, and undergo changes of 
superluminal speed and flux density. The curvature of the tra- 
jectories is strongest within about 2 mas from the core, where 
components C4 and C5 have been traced; C4 has made at least 
one complete oscillation around the mean jet axis. Newly de- 
tected features close to the core also appear to move on differ- 
ent curved paths (Zensus et al. 1995; Krichbaum et al. 1993a). 
At larger distances, the apparent paths merge and straighten 
out. 

These general properties were reconciled by Qian et al. 
(1992) in a geometrical model that explains the trajectories 
of C4 and C5 in 3C 345 as motion along helical paths. They 
introduced smooth (parabolic) bending of the jet axis in the 
sky plane to account for the observed bending of the jet at 
core separations larger than about 3 mas. While this model 
gives a satisfactory fit to the component trajectories and also 
allows derivation of co-moving un-boosted source parameters 



for the origin of the helical structures, and it does not explain 
the straightening of the jet at larger distances from the core. 

In this paper, we discuss the trajectories of components in 
3C 345 by implementing the kinematic aspects of the model of 
Camenzind (Camenzind & Krockenberger 1992). This model 
predicts helical structures as a natural consequence of twisting 
magnetic field lines that are 'frozen' in a wind from a rotating 
accretion disc. From conservation of angular momentum, the 
model also predicts that the curvatures cease at large distances 
from the core. We ignore the complication of overall jet bend- 
ing and restrict ourselves to the inner 3 mas from the core. We 
further assume motion with constant Lorentz factor 7, con- 
stant specific angular momentum i z , and a constant opening 
half-angle ip (in the following called "opening angle") of the 
jet. For generality, we consider other combinations of conser- 
vation quantities in addition to those of Camenzind's model. 
We also include the specific momentum along the jet axis p z , 
which may be used to parameterize other observed jets or mod- 
els (e.g., Hardee 1987). 

Section 2 describes the major variants of our model and the 
analytical solutions to the basic equations of motion. In Section 
3 we apply the model to 3C 345; the results are summarized in 
Section 4. 

2. The Helical Models 

The moving features in the milli-arcsecond jet of the quasar 
3C 345 are modeled as emission regions whose centres of grav- 
ity are moving on 3-dimensionally bent trajectories (Fig. 1). 
The motion is determined by the conservation of the specific 
kinetic energy (equivalent to a constant Lorentz factor 7) of 
the plasma, the specific angular momentum i z with respect 
to the jet axis, the specific momentum component along the 
jet axis p z , and the opening angle of the jet ip (Fig. 2). The 
jet axis is assumed to be a straight line and is oriented in the 
z-direction at an angle 9 to the observer's line of sight. The 
conservation laws for kinetic energy (1), angular momentum 
(2), momentum (3) and jet opening angle (4) can be expressed 
as: 



VLBI-core 
position 




helical path 



Fig. 1. Coordinate system used in the model calculations. 



We now turn to the analytical solutions of the equations of 
motion, which we will use to calculate the position of the jet 
component as a function of time. We start the calculations at 
some initial point (z = 0) at a distance r from the jet axis, 
with some initial angular velocity co and Lorentz factor 70 . 
Initial values are denoted by a subscript 'o'. The origin of the 
coordinate system is identified with the position of the VLBI- 
core of 3C 345 which is assumed to be stationary (Bartel et al. 
1986). 
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Throughout, we use the following notation: 
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Lorentz factor of the plasmon 
opening half-angle of the jet 
inclination angle of the jet axis 
rest mass of the emitting plasmon 
distance from the jet axis 
time in the frame of rest of the quasar 
angular rotation velocity (d(f>/dt) 
cylindrical coordinates (fig. 1) 
velocity in units of the speed of light 
velocity components in units of the 
speed of light 



Fig. 2. Three out of four quantities are assumed to be constant 
during the motion of a jet component: the Lorentz factor 7, the 
specific angular momentum L z and the specific momentum along 
the jet axis p z , and the opening angle of the jet Tp. Three of them 
are needed for every specific case. 



2.1. Case 1 

This model variant assumes the conservation of only mechan- 
ical quantities (momentum along the jet axis p z , the kinetic 
energy Ekin, and the specific angular momentum i z ). There 
are no assumptions for the jet opening angle. Integration of 
the equations of motion (3, 1 & 2, in this order), leads to 



2(t) 
T (t) 

4>( t ) 



V z t + Z 

yjavl t 2 + 8v z t + C 
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arctan 



(5) 
(6) 
(7) 



C = rl, t 



T, = (irl-L*y>*. 



Here and in the following sections the constant of integration 
z is set to z = 0. 

Eqs. (5) & (6) show that the whole jet (or the envelope of all 
possible trajectories) is a rotational hyperboloid of increasing 
or decreasing opening angle (depending on the parameters a, 8, 
and £). It follows from Equation (7) that for this case no com- 
plete revolution of the component around the jet axis is possi- 
ble. Taking into account the properties of the arctan-function 
in the limit t — > oo it follows that <f>oo — 4>o < f • Thus only 
less than a quarter of a revolution can be achieved. The helix 
opens right away, and transforms into an almost straight line 
before undergoing a considerable bend. Therefore, this model 
variant is inadequate for a description of the oscillating motion 
of C4 in 3C 345 (fig. 7). Thus, we conclude that at least one of 
the quantities like the specific angular momentum, the Lorentz 
factor and the momentum component along the jet axis is not 
conserved by the moving emission regions in 3C345. 



2.2. Case 2 

Case 2 is determined by the momentum along the jet axis p z 
the the kinetic energy Ekin and the opening angle ip. 

Integration of the equations of motion (1, 3 & 4) leads to 



2(t) 
T (t) 

4>(t) 



v z t tan^> + r 
<t>° + 



v z tan^> r 



(8) 
(9) 
(10) 



Figure 3 illustrates the general properties of the trajectory, i.e. 
the self-similarity and the constant pitch angle of the helix. 
These properties follow from computing the time derivative 
lo = d<f> I dt of Equation (10). It is found that the azimuthal 
velocity component v$ = tlo is constant in this case. From this 
result and the assumption of a constant momentum along the 
jet axis (3), we also find the pitch angle e of the helix to be 
conserved along the jet: 
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dz 
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const, 



(11) 



This and the assumption of a constant opening angle cause 
an increasing "wavelength" of the helix along the jet (Fig. 
3). This can be seen from calculating the wavelength A( z ) : 
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where k = (v z tan^>)/(w r ). \ z ) an d P(t) are linear in their 
variables and the parameters do not vary along the path. The 
trajectory therefore is self-similar in space and time, i.e. it can 
be scaled linearly in all four coordinates separately. In princi- 
ple this model variant can explain helical motion, but only in 




rel. R.A. [mas] 

Fig. 3. A typical trajectory showing the properties of Case 2. The 
trajectory is self-similar in space and time. 



a range where the helical amplitude is not dampened. 
This case is equivalent to the kinematics of the hydrodynamic 
isothermal helical model presented by Hardee (1987). His set 
of equations r = z^S + r , cf> = 4>o ± 2 ^^° ln(r/r ) (eq. 19a, 
adapted to our notation) and the condition of a constant open- 
ing angle, are similar to our scenario. This can be shown by 
setting 'F = tan ip (or 'F = ip for small opening angles ip) and 
expanding the exponential in Equation (12) into a power se- 
ries. Truncation after two terms yields an initial wavelength 
A ~ 2irVz/io using z = z = 0. Inserting this result into our 
equation for ip^j yields exactly Hardee's Equation 19a. The 
main difference to our simplified model is, that in the hydro- 
dynamical model the jet as a whole is bent. Therefore, it is not 
easy to explain the different paths taken by closely following 
jet components (like C4 and C5), since the overall jet bending 
is assumed to change slowly compared to the motion of the 
individual components. For an application to 3C 345 this case 
is inadequate, since the oscillation of the component path con- 
tinues up to infinite distances. 



2.3. Case 3 

This case is based on the conservation of Ekin, L z , and ip along 
the trajectory and represents an approximation to the kine- 
matics expected from the model by Camenzind (Camenzind 
1986, Camenzind & Krockenberger 1992). 

The solution of the corresponding system of Equations (1,2 
& 4) again provides the position as a function of time: 
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(14) 
(15) 
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where a = vsmtp, b = \J r o — ■jt, c = L/v. Note that ex- 
pression (16) is similar to Equation (7) except for the factor 
1/smip. This counteracts the limit introduced by the arctan 
function on the phase angle 4>(t)- Therefore, in this model the 
number of revolutions is mainly restricted by the opening angle 
if). The limiting opening angle ip allowing at least one complete 
revolution (i.e. Acf> = 2w) about the jet axis — as observed for 
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the second arctan term in Equation (16) may be neglected, 
and for t — > oo we obtain the condition 



A<f> = 2?r = 



sin ^ n 



■ arctan(oo) 



2 sin ^) max ' 



(17) 



It follows that ^/>max — 14.5°. Thus the parameter space for C4 
can be restricted to ip < 14.5°, because at least one oscillation 
has been observed since its discovery (fig. 7). 
We point out that the number of revolutions observed with 
VLBI is not necessarily equal to the number of revolutions 
made by the component since its creation. VLBI-observations 
do not penetrate the jet down to the region where components 
are created. This is supported by optical variability observa- 
tions and theoretical simulations, which imply that most of the 
revolutions take place closer to the central engine (Schramm et 
al. 1993; Camenzind & Krockenberger 1992). In our analysis, 
we start calculations at the position of the VLBI-core of the 
jet. 

Figure 4 shows a "typical" trajectory of this model. In con- 
trast to Case 2 we find no self-similarity, but an asymptotic 
behavior for z — > oo. Close to its origin the coil of the helix is 
strongest, i.e. the pitch angle e is smallest, but it opens along 
the jet axis. These properties match the characteristics of the 
component trajectories observed in 3C345. 



2.4. Case 4 

In this case, the quantities p z , i z and ip are conserved. Again, 
the solutions are derived by integration (now we use Equations 
(3, 4 & 2): 
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For clarity we write 0( t ) as a function of r( t ), which itself is 
a linear function of time (eq. 19). For large z the trajectory 
degenerates into a straight line, since the phase angle 0( t ) ap- 
proaches a fixed value 4>oo, given by 



</>oo = lim <t>( t ) = <t> + 



Lz 



Vz t tan if) 



(21) 



The trajectories of this case are very similar to those in 
Case 3 (Fig. 4). In the following paragraph we discuss the 
discrimination between these two cases and te application to 
3C345. 




rel. R.A. [mas] 

Fig. 4. A typical trajectory illustrating the properties of models 3 
and 4. The trajectory asymptotically transforms into a straight line. 
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Fig. 5. The apparent velocities versus time in the source frame are 
plotted for Cases 3 and 4 corresponding to the path in Fig. 4. The ap- 
parent velocity oscillates strongly where the trajectory shows loops, 
later converging towards a constant value were the helix opens. For 
Case 4 the apparent velocity shows an overall decreasing trend due 
to the non-conservation of the angular momentum and the associ- 
ated loss of the rotational velocity component. In an astrophysical 
jet this could be due to the interaction with the ambient medium. 



This is in accordance with the result that component speeds in 
3C 345 appear to increase with distance from the core (Zensus 
et al. 1995). We note, that the preference for Case 3 here does 
not preclude that the other model variants might be applicable 
in other sources. 



3. Application to 3C 345 

Of the four cases discussed above, we rejected Cases 1 and 2 
as they do not reflect the basic character of the trajectories 
observed in 3C345. Cases 3 and 4 can produce very similar 
trajectories (fig. 4), but they yield different typical velocity 
profiles as shown in Fig. 5. We favor Case 3 over Case 4, as 
it does not only give a better formal fit to the kinematic data, 
but unlike Case 4 it also predicts that component speeds re- 
main on a high level as a component separates from the core. 



The final set of parameters for the chosen model variant was 
derived by fitting to the (different) curved trajectories and ve- 
locity profiles for components C4 and C5 in 3C 345 (Fig. 7-10). 
For our purpose, velocities were computed from subsequent ob- 
servations at a given frequency. Furthermore, we restricted the 
analysis to the inner 3 milliarcseconds from the core. At larger 
distances, the components follow a similar curved path that is 
indicative of large-scale bending of the jet axis not accounted 
for in our model at this point (cf. Unwin & Wehrle 1992, Qian 
et al. 1992, Zensus et al. 1995). 
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Fig. 6. The range covered by the angle to the line of sight determines 
the number of maxima during one turn around the axis on the helical 
path. Ranges 1 and 3 are separated by #max where the apparent 
speed is at maximum. 



3.1. Constraints 

The projected trajectory on the sky plane is fully determined 
by eight partly-interdependent parameters, assuming that the 
origin of the coordinate system is coincident with the sta- 
tionary VLBI-core. They correspond to the three conservation 
quantities kinetic energy, angular momentum, and opening- 
angle; to the three constants of integration; and to two param- 
eters orienting the trajectory on the sky. 

The practical model calculations are based on parameters 
equivalent to the above: the initial value for the Lorentz factor 
7 of a component, the angular velocity co , the opening angle 
if), the time of appearance from the VLBI-core t e] , the distance 
from the jet axis r , the phase angle of the helix <j> , the angle 
between the jet axis and the observer's line of sight 9, and the 
position angle of the jet axis in the sky plane \- in addition, 
we allow for an offset of the observed VLBI-core from the jet 
axis in right ascension (Act) and in declination (AS). 
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Fig. 7. Fits to the observed trajectory of C4 (corresponding to so- 
lutions A, B, and C in Table 1). Filled circles represent the data. 
Where not shown, error bars are smaller than the symbol. 
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Fig. 8. Observed apparent speed of component C4 compared to the 
model calculations (see the trajectories in fig. 7). 



where 7 = (1 — /? 2 ) 2 is the Lorentz factor of the bulk motion, 
and D = [7(1 — /3cos#)] -1 is the Doppler factor. 



Several constraints follow directly from the observations for 
component C4. The lower limit of the Lorentz factor 7 is ap- 
proximately given by the highest observed superluminal speed 
t mal « 5 (fig. 8, Ho = 100 km s" 1 Mpc" 1 , q = 0.5). This im- 
plies that that the angle to the line of sight t9( t ) must be smaller 
than t9( tx ) ftj 20° for these points, if we assume that the Lorentz 
factor is not larger than 100. The opening angle, inferred from 
the envelope of the trajectories, is ipobs ~ 13° (Biretta et al. 
1986). Accounting for projection effects, the true opening angle 
is if) < 13°. The angular speed of rotation around the jet axis 
can be estimated from the observation that C4 has undergone 
a complete oscillation about the axis in t b s ~ 9 years. If 7 fa 5, 
we obtain an estimated lower limit for the initial angular ve- 
locity, corrected for redshift and light travel-time effects (Qian 
et al. 1995) using cos 9 = fi: 
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(22) 



The limit for the initial radial distance from the jet axis r 
follows from the smallest observed core-distance of C4. t e] is 
estimated from the zero separation time as given by Zensus 
et al. (1995). In summary, we restrict the allowed parameter 
space for our model to the following values: 
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The observation that the outer components C3 and C2 fol- 
low almost straight trajectories (Zensus et al. 1995) is used as 
a further constraint. 



Table 1. Parameters of the best fits found for model 3. As indicated in column 2, A and B are solutions of small angles to the 
line of sight (range f in Figure 6) and positive or negative sense of rotation. Solution C is a large-angle solution (range 3 in 
Figure 6) showing different senses of rotation for components C4 and C5. 
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Fig. 9. Fits to the observed trajectory of C5 (see also Table 1). 




Fig. 10. Observed apparent speed of component C5 compared to 
the model calculations (see the trajectories in fig. 9). 



4. Results 

We fitted various sets of parameters and found that only few 
solutions are possible that are qualitatively different. They can 
be distinguished by their sense of rotation (given by the sign of 
Lo ) and the range covered by the angle to the line of sight #( t ) 
during motion along the trajectory (see the (9, « app )-diagram 
in Figure 6). With the three ranges indicated in Figure 6 and 
the two senses of rotation there are six qualitatively different 



solutions possible for every component. Of these, three have 
been found for every component to fit the trajectories (see Ta- 
ble 1). 

In summary we find that : 

(a) The basic motion of components C4 and C5 in 3C 345 
can be explained by a helical path based on the conserva- 
tion of kinetic energy and angular momentum combined 
with a constant opening angle of the jet (Case 3 in Sec- 
tion 2.3). 

(b) Formally the best fit of the individual trajectories and 
the apparent velocities is represented by the large-angle 
solution C (the angle to the line of sight is larger than 
the angle 9 mix of maximum apparent speed !) mix for a 
given Lorentz factor 7). Here components C4 and C5 
have opposite senses of rotation. 

(c) The trajectories of C4 and C5 could be fitted with sim- 
ilar sets of parameters only for small-angle solutions A 
and B (the angle between the jet axis and the line of 
sight was smaller than 9 mix for a given Lorentz factor 
7). Here components C4 and C5 have the same senses of 
rotation (in contrast to the large-angle solution C). We 
regard this as the only reasonable solution considering 
the physical fundamentals of our model. We cannot ex- 
clude, however, that apparent opposite senses of rotation 
might be observable, if the observed apparent motion re- 
flects merely phase velocities or if solid rotation of the 
entire jet is superimposed on the helical motion of indi- 
vidual components. 

(d) The fit is considerably better if the jet axis does not pass 
through the VLBI-core (requiring shifts in both right 
ascension and declination of approximately 0.07 mas). 

(e) The derived Lorentz factor for C5 is always smaller than 
the one obtained for C4 (see Table 1). The deduced dis- 
tance from the jet axis are larger for C5 than for C4. 
Components C6 and C7 seem to move on largely differ- 
ent paths than components C4 and C5 (Krichbaum et 
al., 1993a, Zensus et al., 1995). Such behaviour is also 
consistent with the view of a shell like structure or nested 
conical jets expected from MHD models (e. g., Camen- 
zind, 1986). 

Solution B represents the best fit found for the combined data 
of components C4 and C5. The main difference compared to so- 
lution A showing in the fits to the apparent velocities; around 



ent velocity of C5 is low, whereas the opposite is true around 
1987. However, except for C4 around 1982, the basic behavior 
of solution B agrees with the observation for both components. 
Solution B yields the following physical parameters (see Table 
1): Lorentz factors 7 = 5.8 and 4.6 for C4 and C5, respec- 
tively, and for both components a mean angle to the line of 
sight 9 = 6.8°, using a opening angle ip of 0.95°. These num- 
bers are very close to the values 7 = 7.2±1.0 and 9 = 6.8°±1.5° 
obtained by Unwin & Wehrle (1992) and also to those deduced 
by Zensus et al. (1995). They find a lower limit for the Lorentz 
factor 7 > 8 and the angle to the line of sight 9 < 5.4°. 



4-1. Error estimates 

Fixing all other parameters, the possible ranges for the Lorentz 
factor 7, for the angle to the line of sight 9, and for the year 
of ejection t e] are about 10 percent. The initial phase angle <j> 
of the helix is determined to within 20°, whereas the position 
angle \ can onl y vary by about 2° to yield a fit consistent 
with the data to about la. The initial angular velocity co and 
the radius r are fixed to within less than 0.5degyr -1 and 
0.1 lightyears, respectively. However, the uncertainty is several 
times larger, if one allows co and r to vary simultaneously. 
Another strong dependency exists between 7 and 9. 7 may vary 
by about ±2, and 9 about ±3°. The opening angle ip may vary 
within a few tenths of a degree. For solution C the error ranges 
probably are larger by an estimated factor of 3. In general, the 
error ranges depend on the range covered by the solution in 
the (9, « app )-diagram (Fig. 6). 



5. Conclusion 

Our analysis demonstrates that the basic properties of compo- 
nent trajectories and the apparent velocities for components 
C4 and C5 in 3C 345 can be described by a simple scenario of 
helical motion with a straight jet axis, based on conservation of 
physical quantities. However, there are likely to be additional 
small-scale effects at work which alter the smooth behaviour 
predicted by the model. These effects must be sufficiently large 
to prevent a unique solution that would explain the complete 
kinematic behaviour of even one single component. Similarly 
to the approach by Qian et al. (1992), it was not possible to es- 
tablish a unique sense of rotation about the jet axis. Extension 
of the model to include a curved overall jet axis, and detailed 
monitoring of new components, especially closer than a mil- 
liarcsecond to the core, might eventually allow such a unique 
solution. 

Additional constraints should be obtained from the observed 
flux density evolution of the components at different frequen- 
cies. For example, from the spectral evolution the nature of 
the emission regions can be clarified (Lobanov & Zensus 1994, 
1995), discriminating geometric and kinematic scenarios. We 
still do not know if the moving emission regions are, e. g., tur- 
bulent plasma inhomogeneities or shock waves, which might be 
decided from synchrotron radiation transfer calculations (cf. 
Gomez et al. 1993). Since motion of a specific component is 
typically not completely regular, turbulent motion also might 
play a role causing deviations from the smooth paths predicted 
by our model. 



turbulence-like motion can only account for the smooth overall 
trajectory of individual emission regions. Closer to the core (i. 
e. within less than one milliarcsecond in the case of 3C345), 
such smooth trajectories should prevail, as there the ordered 
magnetic fields associated with the central engine should be 
dominating the source kinematics. 
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